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I.  ZMTRODUCTIOM 


Equatorial ly  trapped  plasmas  are  ion  and  electron 
distributions  trapped  within  a  few  degrees  of  the  Earth's 
magnetic  equator.  Equatorially  trapped  plasmas  can  be 
described  by  a  bi-Maxwellian  distribution  function.  These 
trapped  plasma  distributions  were  defined  by  the  initial 
observations  by  Olsen  (1981) .  Ion  and  electron  distributions 
with  highly  anisotropic  pitch  angle  distributions,  peaked  at 
90  degrees  pitch  angle,  were  observed  at  energies  from  a  few 
eV  to  hundreds  of  eV,  near  geosynchronous  orbit.  These  trapped 
distributions  are  of  interest  as  indications  of  basic  wave 
particle  interactions,  and  as  an  intermediate  process  in 
plasmasphere  filling. 

The  energy  -  pitch  distribution  indicates  the  wave 
particle  interaction  aspect  of  these  plasma  distributions. 
They  are  indicative  of  perpendicular  acceleration  (Tp„p  >  Tp,r)  , 
and  quasi-linear  diffusion  (flat  diffusion  at  low  energy) . 
Though  not  yet  proven,  there  are  indications  of  a 
correspondence  between  equatorially  trapped  plasmas  and 
Bernstein  mode  waves  (equatorial  noise)  and  electron  cyclotron 
harmonics  (Gurnett,  1976  and  Kurth  et  al.,  1979). 

The  plasmasphere  filling  role  is  indicated  by  the 
correspondence  between  the  plasmapause  region  and  the  location 
of  equatorially  trapped  ion  distributions  (Horwitz  et  aJ . , 
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1981) .  The  variations  in  pitch  angle  structure  with  latitude 
also  suggests  the  role  (Olsen  et  al.r  1987). 

There  have  been  previous  surveys  of  equatorial ly  trapped 
plasmas.  Olsen  et  al.,  (1987)  surveyed  DE  1/RIMS  (ion)  data 
for  0  -  100  eV.  Sagawa  et  al.  (1987)  surveyed  the  DE  1/EICS 
(ion)  data  for  0-1  keV.  Both  surveys  of  equatorially  trapped 
plasmas  were  limited  in  altitude  by  the  DE  1  orbit,  which  had 
apogee  at  L  =  4.7.  Both  of  these  surveys  also  lacked 
complementary  electron  data. 

Braccio  (1991)  surveyed  AMPTE/CCE  data  for  both  ions  (30  - 
150  eV)  and  electrons  (150  eV)  out  to  8.8  Earth  radii. 

The  first  purpose  of  this  thesis  is  to  extend  the  survey  done 
by  Braccio  by  conducting  a  statistical  survey  of  the  next  two 
higher  energy  channels  for  both  the  ions  and  electrons  (i.e. 
240  and  442  eV  for  the  ions  and  340  and  770  eV  for  the 
electrons.)  The  survey  will  be  used  to  extend  the  findings  in 
his  paper  and  to  attempt  to  answer  questions  he  raised,  such 
as  whether  the  higher  energy  ions  inhabit  the  region  of 
decreased  probability  from  1200  to  1400  local  time.  The 
electron  survey  will  ensure  that  the  trapped  distribution  is 
not  being  under-evaluated  and  that  the  location  and  shape  of 
the  trapped  distribution  found  by  Braccio  is  accurate. 

The  second  purpose  of  this  thesis  is  to  conduct  a  detailed 
analysis  of  13  individual  days  to  examine  the  nature  of 
equatorially  trapped  distributions.  The  trapped  distributions 
will  be  described  as  bi-Maxwellian  distributions,  in  order  to 
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quantify  a  trapped  distribution.  In  addition,  this  thesis  will 
be  used  to  better  define  the  electron  distributions,  and  the 
location  of  the  trapped  electron  distributions. 
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ZI.  BACXOSODllD 


A.  THE  PLASKASPHEKE 

A  magnetosphere  is  the  region  around  a  magnetized 
planetary  body  in  which  that  body's  magnetic  field  plays  the 
dominant  role  in  defining  the  behavior  of  charged  particles. 
It's  outer  boundary,  the  magnetopause,  occurs  where  the  solar 
wind,  and  the  magnetic  field  in  the  solar  wind,  becomes 
dominant.  This  boundary  occurs  on  the  sunward  side  of  the 
Earth  at  approximately  10  Earth  radii  (roughly  63,750  km) .  The 
location  of  this  boundary  is  determined  by  a  balance  between 
the  pressure  exerted  by  the  solar  wind  and  the  obstacle  formed 
by  the  Earth's  magnetic  field.  During  active  times,  the 
magnetopause  has  been  observed  as  close  as  5  geocentric  Earth 
radii  (Re) .  The  inner  boundary  of  the  magnetosphere  occurs  at 
the  top  of  the  ionosphere.  This  boundary  can  be  taken  as 
occurring  at  an  altitude  of  1000  km  or  1.16  R^.  (Parks,  p.7) 

As  can  be  seen  in  Figure  1,  the  Earth's  magnetosphere  is 
highly  asymmetric.  While  the  sunward  boundary  is  located  at 
approximately  10  Re,  the  Earth's  magnetic  tail  has  been  seen 
to  extend  beyond  200  Re  on  the  nightside.  The  length  and  shape 
of  the  magnetic  tail  again  depends  on  the  interaction  between 
the  geomagnetic  field  and  the  solar  wind.  (Parks,  pp.  7-8) 
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For  our  purposes,  the  major  components  of  the 
magnetosphere  are  the  plasmasphere,  the  plasmasheet,  and  the 
plasmapause.  The  plasmasphere  is  the  region  of  the 
magnetosphere  that  is  closest  to  the  Earth.  It  begins  just 
above  the  ionosphere,  at  low  to  mid-latitudes.  The 
plasmasphere  extends  in  altitude  to  between  3  and  5  in  the 
equatorial  plane,  and  between  +-  60“  magnetic  latitude  just 
above  the  ionosphere.  The  plasmasphere  corotates  with  the 
Earth  and  particles  in  this  region  are  affected  by  the  Earth's 
corotational  Electric  field.  (Parks,  pp.  11  and  73) 

This  region  contains  plasma,  ionized  atoms  and  electrons, 
with  densities  of  10^  -  10^  cm‘^.  Characteristic  ion  and 
electron  energies  are  on  the  order  of  1  eV  at  4.5  Rg.  The 
density  of  the  plasmasphere  decreases  with  altitude.  In 
general,  the  density  in  this  region  experiences  a  gradual  drop 
proportional  to  the  fourth  power  of  the  Mcllwain  L  parameter 
(a  measure  of  altitude  based  on  the  magnetic  field  lines  that 
will  be  discussed  later).  This  is  illustrated  in  Figure  2. 
(Chappell  et  aJ.,  1970) 

At  approximately  3  to  5  Earth  radii,  again  depending  on 
the  magnetic  activity  history,  the  plasmapause  is  encountered. 
This  is  a  transition  region  for  the  plasma  in  which  plasma 
energies  sharply  increase  (Parks,  pp.  231  and  502)  and 
densities  drop,  generally  very  sharply,  that  is  used  to  define 
the  inner  boundary  of  the  plasmapause.  Figure  3. (Harris  et 
al.,  1970). 
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Figure  2.  Plasma  Densi 


Composite  of  the  first  four  H*  ion  concentration  profiles  for  the  inbound  and  out¬ 
bound  passes  of  OGO  5. 


Figure  3.  Plasmapause  Magnetic  Activity  Dependence 
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These  aspects  of  the  plasmasphere  can  be  further 
illu'=trated  using  (relatively)  more  modern  data  from  ISEE  1 
total  electron  density  measurements  obtained  from  observations 
of  plasma  waves.  Figure  4a  shows  density  versus  L,  with  a 
solid  line  at  100  x  (L/4.5)'^  superimposed.  The  plasmapause  is 
at  L  =  4.8,  1700  local  time.  The  plasma  density  outside  the 
plasmapause  continues  to  drop  as  L~*.  This  characteristic  of 
the  plasma  density  profile  can  more  easily  be  seen  if  the  data 
are  normalized  by  V ,  as  in  Figure  4b.  (Olsen,  1992) 

The  plasmasphere  density  is  dependent  on  the  magnetic 
activity.  A  large  magnetic  storm  can  effectively  push  the 
plasmapause  in  to  less  than  3  R^.  Figure  3  shows  the  effects 
of  magnetic  activity  on  the  density  and  location  of  the 
plasmapause.  Magnetic  intensity  increases  from  a  low  in  the 
upper  left  panel  in  the  figure  to  a  maximum  in  the  lower  right 
hand  panel.  (Harris  et  ai.,  1970) 

The  storm-time  electric  field  strips  away  the  plasma  at 
higher  altitudes,  as  the  plasma  is  convected  to  the 
magnetopause.  This  region  is  then  refilled  from  the  ionosphere 
after  the  storm-time  field  relaxes.  The  process,  termed  polar 
wind,  is  driven  by  ambipolar  diffusion  after  the  electric 
field  relaxes  back  to  a  steady  state  value  of  approximately  1 
mV/m.  This  diffusion  process  calls  for  electrons  to  leave  the 
upper  ionosphere,  probably  driven  by  photoemission,  and  move 
along  the  geomagnetic  field  lines.  The  resulting  ambipolar 
electric  field,  caused  by  the  displacement  between  the 
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Figure  4.  Plasma  Density  L 


electrons  and  O*  in  the  upper  atmosphere,  causes  lighter  ions, 
such  as  and  He%  to  be  dragged  up  the  field  lines  after  the 
electrons.  This  'polar  wind*  results  in  a  refilling  rate  of  1 
to  10  ions/cm^  per  da> .  (Horvitz,  1983) 

On  the  nightside  the  plasmasphere  is  bounded  by  the 
plasmasheet,  a  region  of  low  density,  hot  plasma  (with 
densities  on  the  order  of  1  cm'^  and  characteristic  energies 
of  1-10  keV)  .  The  corresponding  plasmapause  for  this  region  is 
very  distinct  and  the  transition  from  plasmapause  to 
plasmasheet  takes  place  rather  rapidly.  This  is  not  the  case 
for  the  region  that  extends  from  just  before  dawn  until  just 
after  dusk  local  time,  on  the  dayside.  (Parks,  pp.  231  and 
502) 

On  the  sunward  side  of  the  Earth,  the  plasmapause  is  a 
region  that  can  be  as  much  as  1  Re  in  width.  Additionally, 
there  is  usually  no  sharp  distinction  corresponding  to  its 
inner  and  outer  boundaries  during  this  local  time  period 
(Parks,  p.231).  Therefore,  it  is  usually  a  matter  of  judgement 
as  to  which  region  you  are  studying. 

The  region  between  the  dayside  plasmapause  and  the 
magnetopause  is  also  ill  defined.  It  is  not  clear  whether  the 
plasmasheet  encircles  the  Earth  and  occupies  this  region. 
While  there  is  no  known  reason  why  this  should  not  be  the 
case,  the  plasma  observed  in  this  region  does  not  display  the 
characteristics  of  that  which  is  found  in  the  nightside 
plasmasheet.  This  has  led  to  questions  concerning  the  plasma 
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filling  mechanism  for  this  region  as  well  as  to  questions  of 
where  the  plasma  in  this  region  comes  from. 

In  the  dusk  region  there  is  an  additional  asymmetry  as 
seen  in  Figure  5.  This  dusk  bulge  is  the  result  of  interaction 
between  the  corotational  electric  field  of  the  plasmasphere 
and  the  cross  tail  electric  field  induced  by  the  solar  wind. 
The  corotational  field  is  the  result  of  the  charged  particles 
rotating  with  the  Earth  while  trapped  in  its  geomagnetic  field 
and  is  directed  radially  inward  toward  the  Earth.  The  cross 
tail  electric  field  is  induced  by  the  solar  wind's  interaction 
with  the  Earth's  magnetic  field.  This  cross  tail  field  is  in 
the  dawn-dusk  direction  in  the  equatorial  plane  of  the 
magnetotail.  The  sum  of  these  two  electric  fields  results  in 
a  series  of  equipotential  contours  which  mirror  the  dusk 
bulge,  Figure  6.  (Parks,  pp.  231-236) 

B.  THEORY 

1.  Plasma  Definition 

A  plasma  is  a  collection  of  discrete  ionized  and 
neutral  particles,  which  has  overall  electrical  neutrality. 
The  physical  dimensions  of  the  plasma  must  be  large  in 
comparison  with  a  characteristic  length  lo  called  the  Debye 
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Equipotentlal  contours  for  the  magnetospheric  electric  field  in  the 
equatorial  plane.  Upper  left:  first-order  approximation  for  die  convection  electric 
field  Be  as  Untform.  The  contours  are  8f»ced  3  kV  apart  hr  Be  =  2.5  x  lO'^V 
m~' .  Upper  right:  the  corotation  electric  field,  contours  spaced  3  kV  apart.  Lower, 
sum  of  convection  and  corotation  electric  fields.  The  heavy  contour  sejrarates  the 
dosed  and  open  convection  regions.  (From  Lyons  and  Williams,  1984) 


Figure  6.  Magnetosphere's  Electric  and  Magnetic  Fields 
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length.  The  total  number  of  charges  in  a  sphere  with  a  radius 
Xd  must  be  much  greater  than  1. 

2 .  Debye  Shielding 

The  electrostatic  potential  of  a  point  charge  g  in  a 
vacuum  is  given  by: 

(1)  — 2 — Volts 

where  f  is  the  electrostatic  potential,  e,,  is  the  permittivity 
constant  and  r  is  the  distance  from  the  point  charge  (Halliday 
and  Resnick,  1988)  .  If  the  charge  is  immersed  in  a  plasma,  a 
positive  charge  will  attract  electrons  while  repelling  ions, 
and  similarly,  a  negative  charge  will  attract  ions  and  repel 
electrons.  The  potential  then  becomes: 

(2)  ,|r=._g_exp'^ 

where  k  is  the  Boltzman  constant,  T,  is  the  electron 
temperature,  which  is  a  measure  of  the  average  kinetic  energy, 
n  is  the  equilibrium  density  of  the  plasma,  and  e  is  the 
charge  of  the  electrons.  This  has  the  effect  of  screening  out 
electric  potentials  in  a  plasma.  The  electron  temperature  is 
used  in  the  definition  of  because  the  electrons  are  more 
mobile  than  the  ions  and  do  most  of  the  shielding  by  creating 
a  surplus  or  deficit  of  negative  charge.  (Parks, 1991) 
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3 .  Plasma  Parameter 


In  order  for  a  collection  of  ionized  particles  to  be 
considered  a  collisionless  plasma,  three  conditions  must  be 
satisfied.  The  Debye  length  must  be  much  less  than  the 
dimensions  of  the  plasma.  The  number  of  particles  in  a  Debye 
sphere  defined  as: 


(3) 


must  be  much  greater  than  1,  that  is,  there  must  be  a  large 
enough  number  of  particles  for  Debye  shielding  to  be 
statistically  valid.  Finally  the  frequency  of  collisions  of 
particles  must  be  low.  The  plasma  parameter  g  is  defined  as: 


(4) 

and  for  a  collisionless  plasma  g  «  1.  (Parks,  1991  and  Chen, 
1983) 

4.  Motion  in  a  Uniform  Magnetic  Field 

A  collisionless  plasma  will  behave  as  a  collection  of 
individual  particles.  The  individual  charged  particles  will 
move  in  trajectories  determined  by  the  applied  electric  and 
magnetic  fields.  For  space  applications  the  fields  resulting 
from  the  particle  motion  are  often  small  and  may  be  neglected. 
This  is  particularly  true  for  the  magnetic  field,  less  so  for 
the  electric  field.  The  force  acting  on  a  charged  particle 
moving  in  a  combined  E  and  B  field  is  given  by  the  Lorentz  equation: 
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(5) 


E=q(E*3SxB) 


where  £  is  the  Lorentz  force,  and  v  is  the  particle  velocity. 

For  the  case  where  £  -  0,  and  the  magnetic  field  is 
uniform,  a  charged  particle  will  execute  simple  cyclotron 
gyration  with  a  frequency: 


(6) 


“c= 


igi-B 

m 


and  radius: 


(7) 


_  ^perp  _  invsintt 
\g\B  ~  \g\B 


where  m  and  v  are  the  mass  and  velocity  of  the  charged 
particle  and  the  pitch  angle  a  is  the  angle  between  the 
velocity  vector  of  the  charged  particle  and  the  magnetic 
field.  The  velocity  parallel  to  the  magnetic  field  line,  Vp„ 
=  vcoso,  is  not  affected  by  the  magnetic  field.  This  motion 
describes  a  circular  orbit  about  a  guiding  center  which  is 
travelling  along  the  magnetic  field  line  with  velocity  Vp„. 
The  trajectory  of  this  particle  is  a  helix  with  it's  axis 
parallel  to  the  field  line.  (Chen,  1983) 

5 .  Magnetic  Mirror 

The  Earth's  magnetic  field  lines  converge  at  both  the 
north  and  south  magnetic  poles,  and  the  field  strength 
increases  with  altitude.  Because  of  this,  there  is  an 
additional  force  that  acts  upon  the  charged  particle.  This 
force  can  be  expressed  as: 
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(8) 


£=-|iVB 


where  m  is  the  magnetic  moment.  Because  the  gradient  in  the 
magnetic  field  is  parallel  to  the  direction  of  the  field  line 
it  can  be  seen  that  this  force  is  directed  along  the  field 
line.  Since  the  force  is  directed  parallel  to  the  particles 
parallel  velocity  component,  it  will  obviously  affect  the 
particle's  velocity.  From  Lentz's  law,  it  can  be  shown  that  n 
is  invariant.  (Parks,  pp.  89-90) 

Since  n  is  invariant,  Vperp  must  increase  as  B  is 
increased.  For  this  to  happen  Vp^^  must  decrease  since  v^perp  = 
v^  -  v^par  (from  conservation  of  energy)  .  Therefore,  given  a 
large  enough  B,  there  will  come  a  point  where  Vpar  ==  0.  At  this 
point  Vp„p  will  equal  v,  and  the  particle  will  mirror  back 
along  the  field  line.  Figure  7.  (Gladstone,  1967).  (Notice 
that  the  gyroradius  also  gets  smaller  as  the  particle 
approaches  the  mirror  point  as  a  result  of  its  Vp^^p 
dependence . ) 

For  a  particle  that  mirrors,  equation  8  then  leads  to: 

2  2 

jgj  ^^perpo  _  ^^^perpm 

QBo  QB„ 

where  subscript  o  refers  to  values  at  the  equator  and  m  to 
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Figure  7.  Path  of  Mirroring  Particle 


those  at  the  mirror  point.  Rearranging  equation  9  gives: 


(10) 


_  ^p9xpo  _  ^pezpo 

v| 


defining  the  pitch  angle  of  a  particle,  a,  to  be  the  angle 
between  velocity  vector  of  the  partide  and  the  magnetic  field 
line  gives  Vp^^p  =  v  sin  a.  Plugging  this  into  equation  10 
gives : 

(11)  :^=sina| 


which  states  that  all  paz^icles  with  a  pitch  angle  Sq  will 
mirror  at  the  location  defined  by  B  =  Bi,.  Particles  with  a  > 
«o  will  mirror  at  lower  altitudes.  (Parks,  pp.  111-112) 

The  magnitude  of  the  Earth's  magnetic  field  is  given 
according  to: 


(12) 


S=S 


08 


V4-3cos^X 

L^cos®A 


where  B^,  is  the  magnitude  of  the  Earth's  magnetic  field  on  the 
Earth's  surface  at  the  magnetic  equator,  X  is  the  magnetic 
latitude,  and  L  is  the  Mcllwain  L  parameter  (Parks,  p.  54). 
The  Mcllwain  L  parameter  is  variable,  given  in  units  of  Earth 
radii,  used  to  label  magnetic  field  lines  with  relation  to 
where  they  cross  the  plane  of  the  magnetic  equator.  Its  value 
is  given  by: 
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(13) 


cos^X 

where  r  is  the  distance  from  the  Earth  center,  in  Re,  to  the 
field  line  at  the  magnetic  equator.  (Par)cs,  p.  115) 

Sxibstituting  equation  12  into  equation  11  gives: 


(14) 


sin^aa= 


cos*la 

^4^3cos^ 


Therefore,  by  defining  an  equatorially  trapped  plasma  to 
mirror  at  a  magnetic  latitude  of  ±  lOo  or  less,  this  requires 
that  a  charged  particle  have  a  pitch  angle  greater  than  69°, 
at  the  equator,  in  order  to  be  equatorially  trapped.  Figure  8. 
It  is  these  trapped  particles  that  will  be  investigated  in 
this  paper. 

6.  Statistical  Distribution 

When  dealing  with  a  system  which  is  composed  of  a  very 
large  number  of  individual  particles  it  becomes  impractical  to 
solve  the  equations  of  motion  for  the  system.  Instead  the 
particles  may  be  treated  statistically  through  the  use  of  the 
distribution  function.  Classically  the  distribution  function 
gives  the  probability  distribution  of  the  values  of  the 
coordinates  and  momenta  of  the  particles.  The  density  of 
particles  in  coordinate  space  is  obtained  by  integrating  the 
distribution  function  over  the  momentum  or  velocity  space. 
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Pitch  Angle  vs  Magnetic  Latitude 


Figure  8.  Pitch  Angle  verses  Mirror  Latitude 
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That  is: 


(15)  ■n=j*.f  (jZ/Jf)  d^v 

where  n  is  the  density,  and  f(£,y)  is  the  distribution  as  a 
function  of  position  and  velocity  integrated  over  three 
dimensional  velocity  space. 

The  two  distribution  functions  of  interest  in  this 
thesis  are  the  bi-Maxwellian,  and  the  Maxwellian  or  isotropic 
distribution  function,  given  by: 


(16) 

(17) 


where  T  p,jp  and  Tpa^.  are  the  characteristic  temperatures  in  the 
perpendicular  and  parallel  directions  with  respect  to  the 
magnetic  field  line. 

C.  PREVIOUS  OBSERVATIONS 

Thermal  plasma  pitch  angle  distributions  seem  to  have  been 
first  studied  by  Horwitz  and  Chappell  (1979)  and  Comfort  and 
Horwitz  (1980)  .  These  authors  used  electrostatic  analyzer  data 
to  study  ion  pitch  angle  distributions  at  geosynchronous 
orbit,  using  ATS  6  data  taken  in  1974.  The  surveys  dealt  with 
data  taken  at  10. 50  off  the  magnetic  equator. 
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Comfort  and  Horwitz  (1980)  observed  two  important  aspects 
of  ion  pancake  distributions  (peak  flux  near  90o  pitch  angle)  . 
The  first  was  the  occurrence  probability  for  the  pancake 
component  of  the  ion  distribution  was  local  time  and  energy 
dependent.  The  highest  probability  of  occurrence  occurred  in 
the  lowest  energy  channel  (20  -  40  eV)  studied  and  for  local 
times  between  1400  and  1800.  Pancake  distributions  wez«?  seen 
42%  of  the  time  in  this  sector  for  ions  of  that  energy. 

The  second  was  that  Comfort  and  Horwitz  observed  that 
field  aligned  ions  and  ions  with  90o  pitch  angle  seem  to  be 
anti-correlated  at  1700.  Figure  9  shows  that  there  is  a 
decrease  in  the  occurrence  probability  of  field  aligned  ions 
when  there  is  a  peak  in  the  pancake  occurrence  probability. 

Horwitz  et  aJ.  (1981)  studied  pancake  distributions  in  low 
energy  (^  100  eV)  ion  data  obtained  from  the  ISEE  1  mass 
spectrometer.  The  H^  distributions  were  often  found  in  the 
vicinity  of  the  plasmapause, Figure  10,  and  usually  just  inside 
the  plasmapause.  Horwitz  et  al.  also  observed  that  the  pancake 
distribution  was  often  seen  in  the  presence  of  colder, 
isotropic  plasma. 

Olsen  (1981)  observed  a  thermal  plasma  population,  trapped 
within  a  few  degrees  of  the  magnetic  equator,  using 
electrostatic  analyzer  data  from  the  SCATHA  satellite.  Figure 
11  shows  the  ion  count  rate,  for  various  ion  energies,  as  a 
function  of  pitch  angle.  The  data  for  this  plot  was  taken  at 
the  equator  at  approximately  1000  local  time  and  5.5  Rg.  This 
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Tottil  percent  occurrence  frequencies 
for  ion  pitch  angle  distributions  having 
deilgnated  components,  as  functions  of  local 
time* 


Figxire  9.  Field  Aligned  and  Pancake  Trapped  Ion  Distributions 
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Ion  pitch  angle  distributions  st  the  e«|ustor  for  dsf  179  of  1979.  Data  are 
plotted  froM  the  FIX  ahd  HI  detectors  at  11,  dt,  193,  523,  and  900  eV.  Fluxes  are 
sealed  by  increasing  factors  to  keep  the*  frote  overlapping.  Maxlaua  values,^ wlth^ 
Increasing  energy  are  1150  c/s,  7300  c/s,  ‘1200  c/s,  550  c/s,  and  5*0  e/s.  The  0* -180* 
range  corresponds  to  looking  sunward,  with  0*  corresponding  to  looking  south. 


Figure  11. 


ion  Pitch  Angle  Distribution 
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figure  clearly  shows  a  trapped  distribution,  centered  at  90© 
and  2700  pitch  angle,  for  ions  of  energies  11  to  103  eV  and, 
to  a  lesser  extent,  for  those  at  523  eV.  The  900  eV  ions  do 
not  show  evidence  of  a  trapped  distribution.  This  figure  also 
shows  a  well  defined  loss  cone  for  the  three  highest  energies. 

Olsen  observed  a  like  distribution  in  the  electron  data. 
Figure  12.  A  source  cone,  centered  at  Qo  and  ISO®  pitch 
angles,  was  seen  in  the  41  eV  electron  flux  concurrent  with 
the  trapped  distribution  at  higher  energies.  This  led  to 
speculation  that  the  field  aligned  particles  were  the 
(ionospheric)  plasma  source,  and  these  particles  were 
subsequently  heated  in  the  traverse  direction.  Note  that  the 
count  rates  in  the  last  two  figures  are  scaled  differently  for 
different  energy  levels  in  order  to  facilitate  presentation  of 
the  data.  Figure  13,  from  Scott  (1991),  shows  a  plot  of  count 
rate  versus  energy  (in  eV) .  The  trapped  electron  distribution 
is  seen  to  exist  in  the  50  to  1000  eV  range,  corresponding  to 
temperatures  of  100  -  200  eV  and  densities  of  1  -  10  cm‘^.  The 
trapped  ions  show  a  peak  in  the  20  to  200  eV  range.  This 
corresponds  to  temperatures  of  20  to  50  eV  and  densities  of  1 

-  10  cm'\  Sagawa  et  al.  (1987)  observed  a  local  time 
dependence  in  the  location  of  the  trapped  ions  in  data  from 
the  Dynamics  Explorer  (DE)  1  satellite.  Sagawa  et  al. 
additionally  saw  that  the  trapped  ions  were  composed  primarily 
of  H*  ions  and  that  these  were  in  the  lowest  energy  bin  (0.01 

-  1  keV)  of  the  DE  1/EICS  summary  plots.  They  reported  that 
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Electron  pitch  angle  distributions  at 
the  equator  for  day  179  of  1979.  Data  are  plot¬ 
ted  from  the  HI  detector  at  41,  523,  and  4730  eV. 
Pitch  angle  conventions  are  as  in  Figure  7. 


Figure  12.  Electron  Pitch  Angle  Distribution 


29 


100.000 


10.000 


3 

o 

u 


•I  I  I  I  I  I  I I  I  — r  I  I  I  I  I'l  IT 

EQUATORIAL  aUXES 
DAY  179  OF  1979 
21:57  TO  21:41  UT 


.  LO  DETECTOR  ELECTRONS 
a-08OTO06<>  , 

CRX2  ! 

f 

1“ 

I" 


m 

•  • 


.1- 


X 


'' 


I 


/ 1'' 

/  I' 


''"''I'l., ''' 

. . 


LO  DETECTOR  IONS 
a-8SOT09SO 


«  FIX  DETECTOR  IONS 
O-90® 


f/  l' 

'll 

,1 


I  \ 

I  I  ^ 

I.  V  . 


'll, 


•l'''!  ,ll 


5*  10  20  SO  100  200  BOO  1.000  2.000 

ENERGY  UV) 


Ion  and  electron  count  rates  as  a  func¬ 
tion  of  energy  from  the  LO  detector  near  90* 
pitch  angle  for  day  179  of  1979.  The  electron 
count  rate  has  been  scaled  by  a  factor  of  2.  The 
count  rates  from  the  FIX  detector  dwells  were  se¬ 
lected  at  their  maxima  (90*  pitch  angle).  The 
difference  between  the  LO  and  FIX  ion  data  re¬ 
flects  degradation  of  the  spiraltrons  for  the  LO 
ion  detector.  The  LO  count  rate  curve  has  been 
traced  and  moved  up  to  overlap  the  FIX  detector 
data  (about  a  factor  of  2).  The  peak  in  ion 
count  rate  at  700  eV  is  a  local  maximum  in  the 
distribution  function  as  well  (see  Figure  10). 


Figure  13.  Energy  Relations  of  Trapped  Plasmas 
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the  Hcllwain  L  value  was  higher,  for  the  peak  ion  occurrence 
probeibility,  in  the  local  noon  and  dusk  sectors  than  it  was 
near  local  midnight,  Figure  14.  Olsen  et  al.  (1987)  also  saw 
this  in  their  statistical  survey.  Olsen  et  al.  noted  that  the 
latitudinal  extent  of  the  high  probability  region  is  local 
time  dependant,  ranging  from  ±  30o  in  the  early  afternoon 
region  to  ±  10°  in  the  early  dawn  region. 

Olsen  et  al.  (1987)  observed,  from  data  collected  by  DE  1, 
that  the  trapped  ion  distribution  was  composed  primarily  of 
H*,  but  that  He*  was  seen  to  have  a  trapped  component,  having 
10%  the  density  of  the  trapped  H*,  approximately  40%  of  the 
time.  In  one  case,  trapped  O*  was  seen  with  a  relative  density 
of  0.1%  that  of  H*.  Additionally,  the  trapped  distribution  was 
observed  to  be  very  localized  about  the  equator.  This  is  seen 
in  the  fact  that  the  ions  change  from  a  field  aligned 
distribution  to  a  trapped  distribution  and  then  back  very 
quickly  as  the  satellite  traverses  the  equatorial  region. 
Figure  15  illustrates  this  aspect  of  the  evolution  in  the 
pitch  angle  distributions. 

Figures  15a,  15b,  and  15c  show  plots  of  flux  verses  pitch 
angle  for  the  magnetic  latitudes  of  -7.9©,  -1.9o 
(approximately),  and  3. 6°  respectively.  In  this  case,  the  He* 
ions  mirror  the  H*  ions,  although  at  about  3.5%  of  its  flux. 
Figures  15d,  15e,  and  15f  show  the  distribution  functions  for 
H*  in  these  time  periods.  Notice  the  drop  in  density  and  the 
increase  in  temperature  as  the  satellite  enters  the  equatorial 
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Occurrence  probability  of  low>enerfy  (0.01-1  keV)  H'*' 
pancake  distribution  with  peak  ion  Bux  above  10^(cm^t  sr)~*, 
within  5^'  of  the  magnetic  equator  for  active  times  (A'p^3-)  as  a 
function  of  MLT  and  L  shell.  L  shell  bin  size  AL=0.5. 


Figure  14.  Trapped  Ion  L  verses  Local  Time  Dependence 
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Figure  15.  Flux  -  Spin  Phase/Density  Fits 
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region.  (Olsen  et  al.,  1987) 

Klumpar  et  al.  (1987)  found  examples  of  eguatorially 
trapped  plasma  in  the  data  from  the  ANPTE/CCE  satellite.  The 
trapped  ion  distribution  was  found  near  the  plasmapause 
interface.  The  temperature  of  these  ions  were  found  to  be  on 
the  order  of  30  -  50  eV.  Like  Olsen  (1981),  they  also  observed 
that  the  angular  distributions  of  the  trapped  ion  distribution 
beccune  narrower  for  increasing  ion  energies. 

Braccio  (1991)  extended  the  survey  of  the  AMPTE/CCE  data 
by  surveying  ions  in  the  30  -  150  eV  energy  channel  and 
electrons  in  the  150  eV  energy  channel  out  to  8.8  Earth  radii. 
In  this  survey  Braccio  concluded  that  the  location  of 
eguatorially  trapped  plasmas  is  species  dependant.  The  ions 
and  electrons  show  a  different  local  time  dependence  in  the 
location  of  their  occurrence  probability  peak.  The  following 
summarizes  his  conclusions: 

Electrons 

•  display  a  uniform  high  probability  distribution  centered 
at  0900  local  time  and  a  L  value  of  6 

•  display  a  weak,  if  any,  L  dependence  on  local  time 

•  trapped  electrons  begin  to  be  seen  at  dawn  leading  to 
speculation  that  their  existence  is  dependent  on 
photoelectron  emission  from  the  Earth's  ionosphere 

•  the  shape  of  the  distribution  is  basically  conical , 
however,  it  drops  off  more  rapidly  than  it  decreases,  with 
respect  to  local  time 
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Ions 


*  distribution  shows  a  strong  L  dependence  on  local  time 

•  high  trapped  ion  probability  region  begins  at  local  dawn, 
for  L  approximately  4,  and  rises  to  a  maximum  at  between 
1400  and  1500  local  time  with  an  L  value  of  8 

•  the  distribution  drops  off  quickly  in  altitude  after 
peaking  as  local  time  increases  there  is  a  region  of 
decreased  probability  in  the  afternoon  sector  for  his  data 
that  he  suggested  may  be  inhabited  by  higher  energy  ions. 
Part  of  this  thesis  will  investigate  this. 

*  the  overall  shape  of  the  high  trapped  ion  prob2d}ility 
region  mirrors  that  of  the  location  of  the  plasmapause 

'  trapped  electrons  seem  to  be  excluded  from  regions  of  high 
trapped  ion  probabilities  and  vice  versa. 

This  thesis  will  further  extend  the  look  at  the  AMPTE  data 
by  conducting  a  statistical  survey  of  ions  in  the  240  eV  and 
442  eV  energy  channels  and  electrons  in  the  340  eV  and  770  eV 
energy  channels. 

There  is  not  always  a  clear  criteria  used  to  define  an 
equatorially  trapped  plasma.  Intuitively,  the  definition  seems 
to  be  that  the  bulk  of  the  plasma  is  confined  to  (mirror 
within)  10°  of  the  magnetic  equator.  One  aspect  of  the 
development  which  follows  is  an  attempt  to  quantify  this 
concept,  and  relate  the  various  definitions,  tests,  and  survey 
techniques.  Braccio  (1991)  defined,  for  the  purpose  of  his 
study,  an  equatorially  trapped  plasma  as  having  a  minimum 
(threshold)  flux  in  the  80®  to  90°  pitch  angle  bin  of  10®  (cm^ 
s  sr)‘^  for  ions  and  5  x  10®  (cm^  s  sr)‘^  for  electrons  with  an 
anisotropy  greater  than  1.5  (The  anisotropy  is  defined  as  the 
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ratio  of  the  fluxes  in  the  80°  to  90o  pitch  angle  bin  to  those 
in  the  60o  to  70o  pitch  angle  bin.)  in  the  same  period.  His 
survey  generally  was  restricted  to  data  taken  within  10°  of 
the  magnetic  equator.  For  the  purpose  of  the  statistical 
survey  in  this  paper,  the  anisotropy  criterion  of  1.5  remains 
the  scune  for  both  ions  and  electrons  and  the  minimiim  flux 
criterion  for  a  equatorially  trapped  plasma  will  be  as 
follows: 

Ions 

•  240  eV:  10‘  (cm^  s  sr)‘^ 

•  442  eV:  10^  (cm^  s  sr)"^ 

Electrons 

•  340  eV  and  770;  5x10*  (cm^  s  sr)’^ 

The  lower  flux  values  for  the  ions  partly  accounts  for  the 
lower  flux  of  ions  at  the  higher  energy  channels. 

Another  criterion  which  could  be  used  to  characterize  a 
trapped  distribution  is  the  ratio  of  the  plasma's 
characteristic  temperatures  in  the  perpendicular  and  parallel 
directions  with  respect  to  the  magnetic  field  line;  Tp^^p  and 
Tpari.  An  equatorially  trapped  distribution  may  be  described  by 
a  bi-Haxwellian  distribution  function  with  these 
characteristic  temperatures  (Scott  1991) .  For  the  purpose  of 
the  case  studies  of  this  thesis  a  ratio  of  Tp^^p  to  Tpa^i  equal 
to  three  will  be  used  as  a  lower  threshold  required  for  an 
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eguatorially  trapped  distribution  to  occur.  The  relationship 
berveen  the  above  two  definitions;  flux  ratios  and  temperature 
ratios,  is  considered  next. 

Applying  the  threshold  of  a  temperature  ratio  of  three 
determines  the  anisotropy  in  flux  that  must  be  met  for  an 
electron  distribution  to  be  considered  trapped.  Figure  16 
shows  plots  of  the  log  of  the  flvix  versus  pitch  angle  for  150 
eV  and  340  eV.  In  each  the  perpendicular  temperature  was  set 
to  75  eV  and  the  parallel  to  25  eV;  a  temperature  ratio  of 
three.  These  temperatures  were  selected  as  representative  of 
values  seen  in  the  data  analysed.  Scott  (1991)  got  slightly 
larger  ratios.  The  density  was  set  arbitrarily  to  1.2  cm'*  as 
a  reasonable  value.  The  density  does  not  affect  the 
illustration. 

In  both  Figure  16a  and  16b  there  are  shown  two  curves.  The 
uppermost  is  the  plot  of  the  bi-Maxwellian  distribution 
function  with  the  above  Tp,tp  and  Tp.^i  values  at  the  magnetic 
equator.  The  lower  of  the  two  curves  is  a  plot  of  the  same 
function  with  the  same  temperatures,  but  mapped  to  a  magnetic 
latitude  of  10'’.  The  peak  flux  of  this  curve  is  equal  to  the 
magnitude  of  the  flux  at  the  magnetic  equator  at  a  pitch  angle 
of  69°.  This  is  because  a  pitch  angle  of  69°  at  the  equator 
maps  to  a  pitch  angle  of  90°  at  a  magnetic  latitude  of  10°. 
With  our  previous  criterion  for  a  trapped  distribution  being 
one  that  is  limited  to  within  10°  of  the  magnetic  equator  this 
establishes  a  lower  limit  on  the  pitch  angles  analyzed  at  the 
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equator  at  69°  and  a  relevant  pitch  angle  range  of  69°  to  90° 
for  the  analysis  of  a  bi-Maxwellian  distribution. 

Taking  the  ratio  of  the  top  curve  to  the  bottom  curve  in 
Figures  16a  and  16b  a  differential  anisotropy  was  derived.  The 
results  of  this  calculation  was  1.67  for  150  eV  and  3.2  for 
340  eV.  The  survey  work  done  here,  and  those  by  Braccio  (1991) 
use  an  integral  anisotropy  calculated  by  taking  the  ratios  of 
the  flux  in  the  80°  to  90°  pitch  angle  bin  to  the  flux  in  the 
60°  to  70°  pitch  angle  bin.  Note  that  the  survey  files  were 
divided  into  pitch  angle  bins  of  10°  in  size.  Using  this 
definition  for  the  anisotropy  we  take  the  ratio  of  the 
integral  of  the  distribution  function  over  pitch  angles  80°  to 
90°  to  the  integral  of  the  distribution  function  over  pitch 
angles  60°  to  70°.  The  results  of  this  is  an  anisotropy  of 
1.953  for  150  eV  and  4.421  for  340  eV.  For  the  purposes  of  the 
detailed  analysis  case  studies  in  this  thesis  the  rounded 
values  of  2.0  and  4.4  will  be  used  for  150  e'^  and  340  eV  , 
respectively,  as  applied  to  a  definition  of  being  bi- 
Maxwellian  and  hence  an  equatorially  trapped  electron 
distribution.  This  will  allow  for  comparison  between  this  work 
and  previous  survey  techniques.  Note  that  the  integral  ratio 
calculated  here  is  higher  than  that  used  in  the  initial 
survey. 
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Figure  16.  Log  of  Flux  verses  Pitch  Angle  -  150  eV  and  340  eV 


Electrons 


D.  THE  EMPTE/CCE  SATELLITE 

The  Active  Magnetospheric  Particle  Tracer  Explorers 
(AHPTE)  mission  consists  of  three  satellites  that  were 
launched  on  August  16,  1984.  The  purpose  of  this  mission 

(Acuna  et  al.,  1985)  was  to: 

1.  investigate  the  transfer  of  mass  and  energy  from  the 
solar  wind  to  the  magnetosphere  and  to  study  its  further 
transport  and  energization  within  the  magnetosphere. 

2.  to  study  the  interaction  between  artificially  injected 
and  natural  space  plasma. 

3.  to  establish  the  elemental  and  charge  composition  and 
dyneunics  of  the  charged  population  in  the  magnetosphere 
over  a  broad  energy  range. 

Two  of  the  satellites,  the  Ion  Release  Module  (IRM)  and 
the  United  Kingdom  Sub-satellite  (UKS) ,  were  concerned 
primarily  with  the  introduction  of  artificially  injected  ions 
into  the  magnetosphere  and  will  not  be  discussed  further.  The 
third  satellite  was  the  Charged  Composition  Explorer  (CCE) , 
Figure  17.  The  purpose  of  this  satellite  was  to  measure  the 
particle  distribution  of  the  naturally  occurring  plasma,  with 
respect  to  species,  energy,  and  pitch  angle,  as  well  as  to 
measure  the  artificially  released  ions  from  IRM.  (Dassoulas  et 
al.,  1985)  The  CCE  was  placed  in  an  elliptical  orbit  around 
the  Earth  with  a  period  of  15.66  hours  and  an  inclination  to 
the  Earth's  equatorial  plane  of  4.8o.  it  had  an  altitude  at 
perigee  of  1108  km  and  at  apogee  of  49,684  km  (roughly  1.2  and 
8.8  R£  respectively)  .  It  was  spin  stabilized  with  a  spin  rate 
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of  10.25  rev/min  (Dassoulas  et  al.,  1985).  Data  began  to  be 
collected  by  this  satellite  on  August  26,  1984. 

The  payload  of  the  CCE,  Figure  18,  consisted  of  five 
experiments;  1}  the  Hot  Plasma  Composition  Experiment  (HPCE) , 
2)  the  Charge  Energy  Mass  Spectrometer  (CHEM)  ,  3)  the  Medium 
Energy  Particle  Analyzer,  4)  the  Magnetometer,  and  5)  the 
Plasma  Wave  Experiment  (Dassoulas  et  al.,  1985).  This  paper 
concerns  itself  with  data  from  the  HPCE  (and,  indirectly,  the 
Magnetometer) . 

E.  THE  HOT  PLASMA  COMPOSITION  EXPERIMENT  (HPCE) 

The  HPCE  consists  of  the  Ion-Mass  Spectrometer  and  the 
Electron  Background  Environment  Monitor  (EBEM) .  The  ion-mass 
spectrometer.  Figure  19,  provides  mass  per  charge  ion- 
composition  measurements  from  very  low  energies  (corresponding 
to  the  spacecraft  potential)  to  approximately  17  keV.  The  ions 
enter  the  detector  through  a  collimator  which  limits  both 
azimuthal  and  elevation  angles  of  acceptance.  The  azimuthal 
limits  are  at  ±5.5®  while  the  elevation  acceptance  angle 
ranges  from  approximately  ±25©  for  ions  at  the  spacecraft 
potential  to  ±7.5©  for  those  at  17  keV.  The  ions  are  sent 
through  a  retarding  potential  analyzer  (RPA)  and  then 
accelerated  through  a  -2960  V  potential.  (Shelley  et  al., 
1985) 


42 


coaimrM  (coi 


>  ICTMOIN  forCNriM.  MMLYZES  IRM) 


CNCWY  mum  (EAt 


umt 


‘■mu.r  SLIT  2  (S2I 


SLIT  S  ISS) 
srEcrai  mi 


T  I  \  ^IMS  AMLTin  (fM) 
J  ^EHEMSr  BETECTM  (CM) 

L  ENER6Y  OETCCTM  (EllI 


ggTtCAL  MTU 


T(y  VIEW 


Figure  19.  The  HPCE  Ion-Mass  Spectrometer 


The  ions  then  pass  through  the  object  slit  and  into  the 
cylindrical  electrostatic  energy  analyzer.  The  electrostatic 
energy  analyzer  is  programmadsle  in  32  energy  per  charge  steps 
from  3  to  20  keV/e.  The  central  portion  of  the  ion  flux  then 
enters  the  mass  analyzer  through  a  second  slit,  with  a  portion 
of  the  spectrum  measured  by  "energy  detectors"  (EDI  and  ED2). 
The  mass  analyzer  consists  of  a  second  cylindrical 
electrostatic  deflection  system  suspended  in  a  978  G  magnetic 
field.  The  ions  that  exit  this  region,  through  the  image  slit, 
are  detected  by  a  high-current  electron  multiplier  (Shelley  et 
ai.,  1985).  This  instnunent  was  active  from  August  26,  1984 
until  it  failed  on  April  4,  1985. 

The  EBEM  consisted  of  eight  independen  18  oo  permanent 
magnet  electron  spectrometers.  Electrons  entered  the  EBEM 
through  a  5^  full  angle  collimator  and  were  then  deflected 
through  180°  by  a  permanent  magnet.  Figure  20.  They  were  then 
focused  onto  an  exit  aperture,  that  defined  the  allowed 
momentum  range,  and  were  then  detected  by  a  channel  electron 
multiplier.  (Shelley  et  aJ.,  1985) 

Both  the  ion  and  electron  data  for  the  AMPTE/CCE  HPCE  were 
processed  into  pool  files.  These  were  the  data  files  used  for 
survey.  These  pool  files  consist  of  data  arranged  in  6.5 
minute  bins  from  0000  to  2400  universal  time.  There  is  a 
separate  data  file  for  each  day's  data  and  each  file  contains 
both  electron  and  ion  data  for  that  day.  (Shelley  et  al., 
1985) 
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For  this  work,  the  ion  flux  measurements  from  the  energy 
detectors  (ED)  are  used,  since  we  were  not  interested  in 
differentiating  between  and  He*  for  this  survey.  The  pool 
data  are  sorted,  by  time,  into  18  logarithmically  spaced 
energy  bins.  The  lowest  four  channels  use  the  "RPA”  mode  data. 
The  bulk  of  the  data  which  are  availcdsle  in  the  pool  file  have 
only  the  fourth  RPA  channel,  which  provides  an  integral 
measurement  from  approximately  30  to  150  eV,  with  a  weighted 
center  at  50  to  65  eV.  The  remaining  channels  extend  up  to 
17keV/e.  The  ion  flux  is  also  sorted  by  time  verses  pitch 
angle,  with  pitch  angle  bins  from  0°  to  90®  in  increments  of 
lOo.  (Shelley  et  aJ.,  1985) 

The  electron  data  is  likewise  sorted  into  8  energy  bins 
from  50  eV  to  25  keV  and  by  pitch  angle  from  0°  to  90®  in  10® 
increments,  each  also  verses  time.  The  energy  channels  for  the 
ion  electrostatic  energy  analyzer  and  for  the  EBEM  are  given 
in  TABLE  I  and  II.  Data  was  also  collected  for  ions  species 
verses  time  verses  energy  and  for  ion  species  verses  time 
verses  pitch  angle  (Shelley  et  aJ.,  1985).  This  aspect  of  the 
data  was  not  used  in  this  paper. 
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TABLE  I.  ION  ENERGY  CHANNELS  IN  THE  HPCE  ON  THE  AMPTE/CCE 


Energy  Channel 

1 

2 

3 

4 


5 

6 

7 

8 

9 

10 
11 
12 

13 

14 

15 

16 

17 

18 


Energy  of  channel 
Center  (keV/e) 

0.0014 

0.0067 

0.020 

0.050 

0.052 

0.055 

0.065 

0.240 

0.442 

0.657 

0.885 

1.127 

1.660 

2.261 

2.941 

3.709 

5.058 

6.668 

9.339 

12.75 

17.11 


Full  Energy  width 
(keV/e) 

0.0028 

0.0076 

0.020 

0.125  RPAl 
0.122  RPA2 
0.115  RPA3 
0.095  RPA4 
0.216 
0.209 
0.222 
0.235 
0.338 
0.567 
0.641 
0.724 
1.056 
1.480 
2.143 
3.041 
3.884 
2.600 


48 


TABLE  II.  ELECTRON  ENERGY  CHANNELS  IN  THE  HPCE  ON  THE 
AMPTE/CCE 


Electron  Detector 

CMEA 

CMEB 

CMEC 

CMED 

CMEE 

CHEF 

CMEG 

CMEH 


Energy  Channel 
Center  (hev/e) 

0.067 

0.150 

0.340 

0.770 

1.74 

3.94 

8.89 

20.1 


Full  Energy  Width 
(kev/e) 

0.051 

0.126 

0.284 

0.643 

I. 45 
3.28 
7.41 

II. 7 


III.  OB8BRVATZOH8 


A.  8TATI8TIC2a.  8URVEY 
1.  Data  Analysis 

Data  collected  from  August  26,  1984  until  December  6, 
1985  were  processed  and  analyzed.  However,  ion  data  were  only 
available  through  April  4,  1985.  This  was  due  to  the  failure 
of  the  ion-mass  spectrometer  on  that  day.  The  stop  date  for 
the  analysis  was  chosen  because  the  satellite  completed  one 
precession  around  the  Earth,  starting  from  August  26,  1984  and 
ending  on  December  6,  1985. 

For  each  day  of  data,  two  different  types  of 
spectrograms  were  produced  for  both  the  ions  and  electron 
data.  The  first  type  was  an  energy  channel-time  spectrogram, 
for  pitch  angles  in  the  80-90  degree  range.  The  second  was  a 
pitch  angle-time  spectrogram.  The  energy  channels  used  were 
the  ion  fifth  and  sixth  channels  (240  and  442  eV)  and  the 
electron  third  and  fourth  channels  (340  and  770  eV) . 

These  spectrograms  were  examined  for  periods  when  the 
satellite  encountered  the  magnetopause.  The  spectrograuns  were 
also  inspected  to  find  periods  of  very  "choppy”  data  and 
periods  when  the  instriiment  was  undergoing  diagnostic  testing. 
Data  in  these  periods  were  then  removed  from  the  data  to 
ensure  that  the  analysis  would  not  be  contaminated  by  such 
artifacts. 
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These  edited  data  files  were  then  surveyed  to  produce 
local  time  vs  Mcllvain  L  probability  distribution  plots.  The 
results  for  ions  and  electrons  were  plotted  separately  then 
compared . 

2.  Local  Time  -  Mcllvain  L  Surveys 
a.  ICN  SURVEY 

Figure  21  shows  a  plot  of  the  probability 
distribution  for  eguatorially  trapped  ions  using  the  lowest 
criteria  from  the  previous  chapter,  namely  the  ion  flux  in  the 
80°“90°  pitch  angle  bin,  at  240  eV,  be  greater  than  10‘ 
ions/(cm^s  sr) ,  that  the  anisotropy  be  greater  than  1.5,  and 
that  the  ions  be  observed  within  10“  latitude  from  the 
magnetic  equator. 

The  grey  scale  for  the  results  runs  from  0%  to  60% 
with  zero  being  white  and  60%  being  black.  The  coverage  plot's 
grey  scale  ranges  from  0  to  200  counts  and  from  white  to  black 
respectively.  The  scales  are  allowed  to  saturate  (peak 
coverage  was  approximately  600  samples  at  apogee) .  Note  that 
the  1800  to  2400  local  time  sector's  zero  occurrence 
probability  is  due  to  lack  of  coverage. 

The  contents  of  Figure  21  are  alternately  presented 
as  a  surface  plot  in  Figure  22.  This  plot  has  x  and  y  axes  of 
local  time  and  L.  The  z  axis  is  probability  of  occurrence,  from 
0  to  100%.  A  contour  plot  is  also  displayed  as  part  of  this 
figure  to  facilitate  the  reading  of  the  surface  plot.  Again, 
the  1800  to  2400  local  time  sector  was  not  sampled. 
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Figure  21.  Trapped  Ions  (240  eV)  -  Flux  gt  10%  Anisotropy  gt 
1.5,  Maglat  gt  10° 


53 


Figure  22.  Trapped  Ions  (240  eV)  -  Flux  gt  10%  Anisotropy  g 
1.5,  Maglat  gt  10°  -  Surface  Plot 


The  high  probability  (greater  than  30%)  region 
these  ions  starts  at  0500  local  time  at  an  L  of  3.5.  As  local 
time  increases  so  does  the  L  value  of  the  peak  probability.  At 
1400  local  time  the  maximxim  L  value,  8,  is  reached.  At  this 
point  the  probability  appears  to  drop  off  sharply  in  L  as 
local  time  moves  toward  dusk.  This,  however,  may  be  an 
artifact  brought  about  by  the  low  coverage  after  1800  local 
time. 

In  Figures  23  and  24  the  results  of  raising  the 
selection  criteria  for  a  trapped  ion  distribution  are  shown. 
In  Figure  23  the  selection  criteria  are  flux  greater  than  10% 
anisotropy  greater  than  1.5,  and  measurements  within  5°  of  the 
magnetic  equator.  In  Figure  24  the  selection  criteria  are  flux 
greater  than  5  x  10%  anisotropy  greater  than  2.0,  and 
measurements  within  5®  of  the  magnetic  ernator.  The  local  time 
dependence  is  similar  in  these  three  figures. 

Figure  25  shows  a  plot  of  the  probability 
distribution  for  equatorially  trapped  ions  for  the  sixth 
energy  channel  (i.e.  442  eV) .  Here  the  selection  criteria  is 
a  flux  greater  than  10*  ions/ (cm*  s  sr) ,  an  anisotropy  greater 
than  1.5,  and  a  magnetic  latitude  between  within  5“  of  the 
magnetic  equator.  The  grey  scales  are  the  same  as  in  the 
previous  two  figures.  Figure  26  is  a  surface  plot  of  the  same 
data. 

The  high  probability  (greater  than  30%)  region  for 
these  ions  starts  at  0500  local  time  at  an  L  of  3.5.  Again,  as 
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Figure  23.  Trapped  Ions  (240  eV)  -  Flux  gt  10%  Anisotropy  gt 
1.5,  Maglat  gt  5° 
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Figure  24.  Trapped  Ions  (240  eV)  -  Flux  gt  5x10%  Anisotropy 
gt  2.0,  Maglat  gt  5° 
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Figure  25.  Trapped  Ions  (442  eV)  -  Flux  gt  10%  Anisotropy  gt 
1.5,  Maglat  gt  5° 


58 


as  local  tiime  increases  so  does  the  L  value  of  the  peak 
prob2d>ility.  At  1400  local  time  the  maximum  L  value,  8,  is 
reached  and  then  the  probability  drops  off  shaxply,  probedsly 
due  to  the  lack  of  data  after  1800. 

The  survey  of  the  additional  energy  channels  for 
the  ions  did  not  dramatically  change  the  results  found  by 
Braccio.  Figures  27  and  28  are  the  probability  distribution 
plot  and  surface  plot  from  Braccio  (1991)  for  ions  in  the 
fourth  RPA  energy  channel  (30  -  150  eV)  with  the  criteria  of 
flux  in  the  80°-90°  pitch  angle  bin  be  greater  than  10® 
ions/(cm^s  sr) ,  that  the  anisotropy  be  greater  than  1.5,  and 
that  the  measurements  be  within  10°  latitude  of  the  magnetic 
equator.  It  is  found  that  the  high  probability  region  (greater 
than  45%)  for  the  ions  started  at  0500  local  time  and  in  the 
area  of  L  =  3.5.  His  maximtim  L  value  is  on  average  about  8.0 
and  occurs  generally  around  1400  local  time.  At  the  maximum 
the  probability  appears  to  drop  off  sharply  in  L  as  local  time 
moves  toward  dusk,  similar  to  the  results  above. 

The  apparent  decrease  in  probability  of  occurrence 
for  ions  in  the  region  of  1200  to  1400  local  time  found  by 
Braccio  is  readily  evident  in  both  energy  channels  surveyed, 
which  seems  to  dampen  Braccio 's  hypothesis  that  the  higher 
energy  ions  inhabited  this  region. 
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Figure  27.  Trapped  Ions  (30-150  eV)  -  Flux  gt  10^,  Anisotropy 
gt  1.5,  Maglat  gt  10® 
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Figure  28.  Trapped  Ions  (30-150  eV)  -  Flux  gt  10*^,  Anisotropy 
gt  1.5,  Maglat  gt  10®  -  Surface  Plot 
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b.  ELECTRON  SURVEY 


Figure  29  shows  the  probability  distribution  for 
trapped  electrons  with  an  energy  of  340  eV  meeting  the 
criteria  of  flux  greater  than  5  x  10*  electrons/ (cm*  s  sr) , 
anisotropy  greater  than  1.5,  and  measurements  within  10°  of 
the  magnetic  equator.  Figure  30  shows  the  seune  type  of  plot 
for  770  eV  electrons  meeting  the  same  criteria.  Figure  31 
presents  this  770  eV  data  as  a  surface  plot.  It  is  readily 
apparent  that  the  electron  distribution  is  vastly  different 
from  that  of  the  ions.  There  is  no  obvious  L  versus  local  time 
trend  for  the  electron  probability  distribution,  similar  to 
the  results  found  by  Braccio,  who  also  fowmd  the  electron 
probability  distribution  to  be  conical  in  shape  with  the 
region  of  highest  probability  being  between  1000  and  1100 
local  time  and  for  L  values  around  6.5  -  7.0,  as  shown  in  the 
probability  and  surface  plots.  Figures  32  and  33,  taken  from 
Braccio  (1991) .  Here,  the  criteria  is  flux  greater  than  5  x 
10*  electrons/ ( cm*  s  sr)  ,  anisotropy  greater  than  1.5,  and 
measurements  within  10°  of  the  magnetic  equator. 

The  cone  is  not  completely  symmetrical.  The 
probability  of  occurrence  increases  to  its  peak  value  more 
gradually,  as  a  function  of  local  time,  than  it  decreases. 
This  characteristic  does  not  seem  to  alter  even  when  changing 
the  selection  criteria  or  energy  channels  and  is  the  same  as 
that  found  by  Braccio.  The  statistical  analysis  of  the 
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Figure  29.  Trapped  Electrons  (340  eV)  -  Flux  gt  5x10*, 
Anisotropy  gt  1.5,  Maglat  gt  10° 
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Figure  30.  Trapped  Electrons  (770  eV) 
Anisotropy  gt  1.5,  Maglat  gt  10° 
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Figure  32.  Trapped  Electrons  (150  eV)  -  Flux  gt  5x10' 
Anisotropy  gt  1.5,  Maglat  gt  10® 
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Figure  33.  Trapped  Electrons  (150  eV)  -  Flux  gt  5x10 
Anisotropy  gt  1.5,  Maglat  gt  10®  -  Surface  Plot 
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coarsely  gridded  survey  files  is  now  extended  using  more 
detailed  looks  at  the  raw  data  files. 

B.  DETAILED  ANALYSIS  CASE  STUDIES 
1.  Data  Analysis 

Detailed  analysis  of  the  electron  data  was  conducted, 
using  data  at  150  eV  and  340  eV.  ion  data  were  not  readily 
available.  Data  were  collected  for  the  days  in  TABLE  III.  Only 
days  84/336  and  85/002  will  be  presented  in  the  thesis.  The 
days  listed  in  TABLE  III  were  selected  from  the  set  of  days 
where  ••large”  electron  events  occurred,  as  found  in  Braccio 

Table  III.  DETAILED  ANALYSIS  DATA  COLLECTION  DAYS 


Day 

Local  Time 
Range 

UKp 

D.t 

84/266 

0941-1634 

26 

-21 

84/283 

0510-1538 

26 

-24 

84/296 

0755-1420 

39 

-55 

84/315 

0655-1349 

23 

-11 

84/336 

0533-1249 

22 

-24 

85/002 

0359-1418 

24 

-3.2 

85/039 

0129-0834 

28 

-24 

85/227 

1440-1938 

20 

-21 

85/238 

1330-2330 

22 

-10 

85/264 

2002-2156 

29 

-19 

85/301 

1618-2049 

7 

-17 

85/320 

1116-1837 

19 

-18 

85/337 

1021-1742 

17 

-31 
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(1991)  and  extended  here,  to  obtain  reasonably  complete  local 
time  coverage.  Special  emphasis  was  given  to  the  morning 
sector  (0600  to  1200  local  time)  .  Each  segment  of  data 
considered  covered  the  period  between  plasmapause  encoiinters, 
roughly  15  hours  of  data  per  day.  Also,  listed  with  the  days 
and  times  in  TABLE  III  are  the  magnetic  activity  indices,  EKp 
and  Djt,  for  the  respective  time  period. 

Data  were  initially  considered  in  one  hour  increments. 
For  each  hour  and  each  of  the  energy  channels  surveyed,  150 
and  340  eV,  a  pitch  angle  vs  UT  time  spectrogram  was 
generated.  Figure  34  is  a  compilation  of  the  one  hour  pitch 
angle  vs  time  spectrograms  in  reduced  form  for  150  eV 
electrons  on  day  84/336.  The  150  eV  spectrograms  for  the 
additional  days  surveyed  can  be  found  in  the  Appendix.  The 
grey  scale  for  the  spectrogram  runs  from  7.2  to  9.0  on  a  log 
scale  of  the  flux.  Pitch  angle  runs  from  0'  -  180°  and  UT  from 
0000  to  1400.  Also,  plotted  on  the  x-axis  is  the  corresponding 
local  time  (LT) ,  Mcllwain  L  value  and  magnetic  latitude.  Note 
the  high  flux  at  90  degrees  from  1000  to  1400  UT.  Note  how  the 
equatorially  trapped  distribution  fades  as  the  satellite 
approaches  perigee  (the  plasmapause  at  L  "  3.9.)  This 
characteristic  can  be  seen  in  a  number  of  the  survey 
spectrograms  as  shown  in  Appendix.  The  data  were  analyzed  by 
considering  the  log  of  distribution  function  as  a  function  of 
cos^o.  Figure  35.  This  functional  form  should  produce  straight 
line  segments,  if  the  data  are  well  represented  as  a 
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Figure  34.  Pitch  Angle  verses  Time  Spectrogram  -  Day  84/336 
Energy  150  eV 
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Figure  35.  Electron  Pitch  Angle  Distribution 
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bi-Maxwellian  (Scott,  1991.)  From  these  segments  values  for 
Tperp  and  Tpari  Can  be  estimated,  and  various  fliix  ratios 
obtained.  Much  of  this  work  is  based  on  consideration  of  how 
the  data  (log  f)  vary  with  pitch  angle  (cos^a)  . 

From  this  figure  a  distinct  break  in  the  vicinity  of 
60  degrees  can  be  seen.  At  lower  pitch  angles  the  isotropic 
background  dominates  the  pitch  angle  distribution.  Because  of 
this  effect  the  data  analysis  range  was  limited  to  between  60° 
and  90°,  a  restriction  that  more  accurately  considers  the 
phenomenon  of  interest.  This  restriction,  however,  has  the 
effect  of  limiting  the  data  set  which  can  be  analyzed.  These 
equatorial  pitch  angles  (60°  to  90°)  do  not  map  far  from  the 
equator.  A  restriction  of  +-6°  magnetic  latitude  was  chosen, 
which  allows  a  suitable  range  of  (equatorial)  pitch  angles, 
and  allows  a  reasonable  range  of  events. 

Fluxes  were  obtained  in  three  ways; 

1.  Average  Flux  (60°  -  70°) :  the  flux  in  the  60°  to  70°  pitch 
angle  bin  at  the  location  of  the  satellite. 

2.  Average  Flux  (80°  -  90°)  ;  the  flux  in  the  80°  to  90°  pitch 
angle  bin  at  the  location  of  the  satellite. 

3.  Fitted  Flux  (60°  -  90°);  the  fitted  flux  was  derived  by 
taking  a  least  squares  fit  of  the  data  at  pitch  angles 
between  60°  and  90°  on  a  plot  of  the  log  of  the  distribution 
function  versus  cos^a.  This  fit  was  used  to  project  the 
values  which  would  have  been  obtained  at  the  magnetic 
equator,  effectively  assuming  a  bi-Maxwellian  form. 
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Two  different  anisotropies  were  calculated  and  are 
defined  as  follows: 

1.  Average  Anisotropy:  the  ratio  of  the  average  fluxes  in 
the  80°  to  90°  pitch  angle  bin  to  those  in  the  60°  to  70° 
pitch  angle  bin. 

2.  Equator  Anisotropy:  the  fluxes  at  the  latitude  of 
measurement  mapped  to  the  equator  and  a  ratio  taken  for  the 
flvix  at  90°  versus  the  flux  at  69°. 


Two  temperatures  were  calculated  for  each  energy 

channel : 

i*  *^perp  •  the  characteristic  temperature  in  the 
perpendicular  direction,  with  respect  to  the  magnetic  field 
line,  of  a  bi-Maxwellian  distribution.  Calculated  from  the 
difference  in  the  log  of  the  distribution  functions  for  150 
eV  and  340  eV  at  a  pitch  angle  of  0°.  The  calculated  value 
is  the  same  for  both  energy  channels. 

2.  Tpgji  :  the  characteristic  temperature  in  the  parallel 
direction,  with  respect  to  the  magnetic  field  line,  of  a 
bi-Maxwell ian  distribution.  This  is  based  on  the  calculated 
Tperp  2ind  the  slope  of  the  pitch  angle  distribution  at  the 
respective  energy.  The  calculated  value  is  different  for 
each  energy  channel,  reflecting  the  non-Maxwell ian  element 
of  the  energy  distribution. 


Both  temperatures  and  the  density  were  calculated  from 
a  bi-Maxwellian  form  for  the  distribution  function.  This 
produces  valid  and  well  behaved  results  if  the  distribution 
can  be  characterized  as  a  bi-Maxwellian.  If  it  is  not  well 
described  as  a  bi-Maxwellian  the  results  will  be  invalid  and 
not  well  behaved. 

The  ratio  of  Tp^rp  to  Tpa^-j  was  also  calculated  for  the 
data.  This  ratio  was  used  to  further  characterize  a  trapped 
distribution.  For  the  purposes  of  this  paper  the  ratio  had  to 


exceed  three  for  a  distribution  to  be  characterized  as 
eguatorially  trapped.  This  corresponds  to  a  density  decrease 
factor  of  20%  at  10°  magnetic  latitude  in  the  aibsence  of 
parallel  electric  fields.  (Scott,  1991) 

2.  Day  84/336  Case  Study  Trapped  Distributions 
a.  150  eV 

Figure  36  shows  the  plots  of  the  previously  defined 
fluxes  vs  UT  time  at  150  eV.  The  solid  horizontal  line  plotted 
with  the  fluxes  in  the  bottom  two  plots  is  for  the  flux  value 
of  5x10*  (cm^  s  sr)*^  ,  the  minimum  threshold  established  in  the 
last  chapter.  Also,  plotted  on  the  x-axis  are  the  local  time, 
Hcllwain  L  and  the  magnetic  latitude. 

Figure  37  shows  the  plots  of  the  previously  defined 
anisotropies  versus  time  at  150  eV.  The  solid  horizontal  line 
plotted  on  each  is  the  anisotropy  threshold  of  2.0,  defined  in 
the  previous  chapter  as  the  minimum  anisotropy  allowed  to 
characterize  a  trapped  distribution.  Note  the  periods  that  the 
anisotropies  exceed  the  threshold  correspond  to  the  periods  of 
increased  flux.  Also,  shown  in  Figure  37  is  the  plots  of  Tperp 
vs  time,  Tpa^i  vs  time,  and  the  ratio  Tperp  to  Tpan  vs  time.  The 
solid  horizontal  line  plotted  on  the  plots  for  temperature 
ratios  vs  time  is  the  ratio  threshold  level  of  three  defined 
in  the  last  chapter.  At  the  time  the  temperature  ratio  exceeds 
three  perpendicular  temperature  is  on  average  70  eV  while  the 
parallel  temperature  drops  from  values  of  20  eV  to  10  eV.  The 
peak  of  the  ratio  has  a  value  of  6.0.  For  this  day  a 
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Figure  36.  Fluxes  verses  Time  -  Day  84/336,  Energy  150  eV 
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eguatorially  trapped  electron  distribution  appears  to  occur 
for  a  Mcllwain  L  value  of  about  6  to  7.5.  Note  that  the  lower 
anisotropy,  1.5  (the  dashed  line  in  Figure  37),  used  in  the 
statistical  survey  would  give  a  wider  L  extent,  from  L  =  3  to 
9.  With  this  lower  anisotropy  and  the  corresponding  L  range  as 
a  boundary,  for  where  eguatorially  trapped  electron 
distributions  exist,  a  flux  value  of  5x10^  and  a  temperature 
ratio  of  two  can  be  interpreted  as  new  criterion  thresholds 
for  eguatorially  trapped  electrons,  as  seen  by  the  dashed 
lines  in  the  flux  and  temperature  ratio  plots  of  Figures  36 
and  37. 

Figure  38  shows  the  plot  of  the  average  anisotropy 
versus  the  eguator  anisotropy  and  the  plot  of  temperature 
ratio  versus  eguator  anisotropy.  The  average  anisotropy  is 
slightly  higher  than  the  fitted  eguatorial  anisotropy  as  can 
be  seen  when  compared  to  a  line  of  slope  one  (the  line  on  the 
plot.)  In  the  plot  of  temperature  ratio  vs  eguator  anisotropy 
the  data  starts  to  exceed  the  ratio  criterion  threshold  of 
three  at  an  eguator  anisotropy  of  about  2.0,  if  the  high  end 
of  the  data  points  is  taken.  This  tends  to  support  our 
definition  for  an  eguatorially  trapped  electron  distribution 
using  flux  ratios,  and  shows  that  a  bi-Maxwellian  description 
gives  comparable  results. 

Using  the  criteria  for  eguatorially  trapped 
electron  distributions  defined  in  the  previous  chapter  of  a 
minimum  flux  for  150  eV  electrons  of  5x10®  (cm^  s  sr)S  an 
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Figure  38.  Average  Anisotropy  and  Temperature  Ratio  verses 
Equator  Anisotropy  -  Day  84/336,  Energy  150  eV 


average  anisotropy  of  2.0,  a  measurement  within  10°  of  the 
magnetic  equator,  and  the  additional  criterion  for  a  ratio  of 
Tperp  vs  Tp„i  to  be  greater  than  three  an  equatorially  trapped 
plasma  can  be  seen  to  start  at  around  0942  local  time  at  a  L 
value  of  8.7  and  at  magnetic  latitude  of  -2.4°.  As  L  decrease 
with  time  all  anisotropies  increase  to  a  peak  at  around  1050 
local  time  and  a  L  value  of  6.9  and  a  temperature  ratio  of 
6.0.  All  anisotropies  decrease  with  time  after  that  until  the 
trapped  distribution  ends  at  around  1147,  L  of  5.30  and 
magnetic  latitude  of  2.4°. 
b.  340  eV 

Figure  39  shows  the  plots  of  the  previously  defined 
fluxes  vs  UT  time  at  340  eV.  The  solid  horizontal  line  plotted 
with  the  fluxes  in  the  bottom  two  plots  is  for  the  flux  value 
of  5x10*  (cm^  s  sr)'^  ,  the  minimum  threshold  established  in  the 
previous  chapter.  Also,  plotted  on  the  x-axis  are  the  local 
time,  Mcllwain  L  and  the  magnetic  latitude. 

Figure  40  shows  the  plots  of  the  previously  defined 
anisotropies  versus  time  at  340  eV.  The  solid  horizontal  line 
plotted  on  each  is  the  anisotropy  threshold  of  4.4,  defined  in 
the  previous  chapter  as  the  minimum  anisotropy  allowed  to 
characterize  a  trapped  distribution.  Also,  shown  in  Figure  40 
is  the  plots  of  Tpe,,p  vs  time,  Tpan  vs  time,  and  the  ratio  Tpe^p 
to  Tpari  vs  time.  The  solid  horizontal  line  plotted  on  the 
plots  for  temperature  ratios  vs  time  is  the  ratio  threshold 
level  of  three  defined  in  the  last  chapter.  At  the  time  the 
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Figure  39.  Fluxes  verses  Time  -  Day  84/336,  Energy  340  eV 
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temperature  ratio  exceeds  three  perpendicular  temperature  is 
on  average  70  eV  while  the  parallel  temperature  drops  to 
values  of  23  eV.  At  this  energy  a  trapped  electron 
distribution  appears  to  occur  for  a  more  limited  Mcllwain  L 
range,  at  the  higher  end  of  the  L  range  defined  by  the  150  eV 
data. 

Using  the  criteria  for  trapped  electron 
distributions  defined  in  the  previous  chapter  of  a  minimum 
flux  for  340  eV  electrons  of  5x10®  (cm^  s  sr)'S  an  average 
anisotropy  of  4.4,  a  measurement  within  10°  of  the  magnetic 
equator,  and  the  additional  criterion  for  a  ratio  of  Tp^^p  vs 
Tpari  to  be  greater  than  three  an  equatorial ly  trapped  plasma 
can  be  seen  to  start  at  around  1031  local  time  at  a  L  value  of 
7.5  and  at  a  magnetic  latitude  of  0.8°.  The  trapped 
distribution  does  not  extend  very  far  ending  at  1041  at  L 
equal  to  7.2  and  magnetic  latitude  of  1.2°.  A  more  realistic 
anisotropy  ratio  of  ’  2  and  a  temperature  ratio  of  2  ,  as 
shown  by  the  dashed  line  in  both  plots  of  Figure  40  would  give 
a  latitude  extent  more  in  line  with  respect  to  the  150  eV 
results . 

3.  Day  85/002  Case  study  Trapped  Distributions 
a.  150  eV 

Figure  41  shows  the  plots  of  the  150  eV  fluxes  vs 
time  for  85/002.  The  solid  horizontal  line  plotted  with  the 
fluxes  in  the  bottom  two  plots  is  for  the  flux  value  of  5x10® 
(cm^  s  sr)''.  Also,  plotted  on  the  x-axis  are  the  local  time. 
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Mcllwain  L  and  the  magnetic  latitude.  Comparison  with  the 
anisotropy  suggests  a  flux  level  of  5x10^  would  be  a  more 
accurate  indicator  on  this  day,  as  shown  by  the  dashed  line  in 
the  plot. 

Figure  42  shows  the  plots  of  the  previously  defined 
anisotropies  versus  time  at  150  eV.  The  solid  horizontal  line 
plotted  on  each  is  the  anisotropy  threshold  of  2.0,  defined  in 
the  previous  chapter  as  the  minimum  anisotropy  allowed  to 
characterize  an  eguatorially  trapped  distribution.  Note  the 
periods  that  the  anisotropies  exceed  the  threshold 
corresponding  to  the  periods  of  increased  flux  (>  5x10’)  . 
Also,  shown  in  Figure  42  are  the  plots  of  Tperp  vs  time,  Tpa^  vs 
time,  and  the  ratio  Tperp  to  Tpari  vs  time.  The  solid  horizontal 
line  plotted  on  the  plots  for  temperature  ratios  vs  time  is 
the  ratio  threshold  level  of  three.  At  the  time  the 
temperature  ratio  exceeds  three  perpendicular  temperature  is 
on  average  80  eV  while  the  parallel  temperature  drops  to 
values  of  20  eV.  For  this  day  a  trapped  electron  distribution 
appears  to  occur  for  a  Mcllwain  L  value  of  about  5  to  7.  Note 
the  lower  anisotropy  ratio  of  1.5  (the  dashed  line  in  Figure 
42) ,  used  in  the  statistical  survey  would  give  a  wider  L 
extent,  from  L  =  5  to  8.  With  this  lower  anisotropy  and  the 
corresponding  L  range  as  a  boundary,  for  where  eguatorially 
trapped  electron  distributions  exist,  a  flux  value  of  5x10’ 
and  a  temperature  ratio  of  two  can  be  interpreted  as  new 
criterion  thresholds  for  eguatorially  trapped  electrons. 
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Figure  41.  Fluxes  verses  Time  -  Day  85/002,  Energy  150  eV 
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as  seen  by  the  dashed  lines  in  the  flux  and  temperature  ratio 
plots  of  Figures  41  and  42. 

Using  the  criteria  for  trapped  electron 
distributions  defined  in  the  previous  chapter  of  a  minimum 
flux  for  150  eV  electrons  of  5xl0‘  (cm^  s  sr)'‘,  an  average 
anisotropy  of  2.0,  a  measurement  within  10®  of  the  magnetic 
equator,  and  the  additional  criterion  for  a  ratio  of  Tp^rp  vs 
Tpari  to  be  greater  than  three  a  equatorially  trapped  plasma 
can  be  seen  to  start  at  around  0900  local  time  at  a  L  value  of 
7.2  and  at  magnetic  latitude  of  2.6°.  As  L  decreases  with 
time  the  anisotropies  increase  to  a  peak  at  around  0909  local 
time  and  an  L  value  of  6.9.  The  anisotropies  decrease  with 
time  after  that  until  the  trapped  distribution  ends  at  around 
1003,  L  of  5.30  and  magnetic  latitude  of  3.6“. 
b.  340  eV 

Figure  43  shows  the  plots  of  the  fluxes  vs  time  at 
340  eV.  The  solid  horizontal  line  plotted  with  the  flxixes  in 
the  bottom  two  plots  is  for  the  flux  value  of  5x10®  (cm*  s 
sr)*.  Also,  plotted  on  the  x-axis  are  the  local  time,  Mcllwain 
L  and  the  magnetic  latitude. 

Figure  44  shows  the  plots  of  the  anisotropies 
versus  time  at  340  eV.  The  solid  horizontal  line  plotted  on 
each  is  the  anisotropy  threshold  of  4.4.  Also,  shown  in 
Figure  44  are  the  plots  of  Tp„p  vs  time,  Tp^^i  vs  time,  and 
theratio  Tp^rp  to  Tpa^i  vs  time.  For  the  short  interval  where  the 
temperature  ratio  exceeds  three  the  perpendicular  temperature 
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Figure  43.  Fluxes  verses  Time  -  Day  85/002,  Energy  340  eV 
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is  on  average  90  eV  while  the  parallel  temperature  drops  to 
values  of  30  eV.  Again,  the  eguatorially  trapped  distribution 
is  much  more  limited  in  its  appearance  at  340  eV.  Indeed  the 
criterion  of  a  ratio  of  4.4  is  not  met.  However,  if  we  lower 
the  anisotropy  criterion  to  3.0  (the  dashed  line  in  the 
anisotropy  plots  of  Figure  44)  a  trapped  distribution  would  be 
defined  to  occur  around  a  Mcllwain  L  value  range  of  6  to  7.5 
at  around  0901  local  time  at  a  L  value  of  7.1  and  at  magnetic 
latitude  of  2.8°.  The  trapped  distribution  does  not  extend 
very  far  ending  at  0928  at  L  equal  to  6.3  and  magnetic 
latitude  of  3.4°.  This  corresponds  to  a  temperature  ratio  '  2 
(the  dashed  line  in  the  temperature  ratio  plot  of  Figure  44) 
4.  Detailed  Analysis  Trapped  versus  Non-Trapped 

Distributions 

In  Figure  34  there  was  evidence  of  two  areas  of 
trapped  electron  distributions  for  day  84/336.  The  first  was 
presented  using  Figures  36  through  38,  which  was  an 
eguatorially  trapped  distribution  as  defined  by  our  criteria. 
The  other  area  that  presents  itself  in  Figure  34  is  the  time 
period  between  0553  and  0700  local  time.  The  broader 
distribution  seen  earlier  in  the  day  fail  our  criteria.  Figure 
45  and  46  show  plots  like  previously  of  the  fluxes, 
anisotropy,  temperatures  and  temperature  ratios,  all  verses 
time  for  the  entire  orbit  on  84/336  at  150  eV.  Though  both 
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time  periods  show  similar  increases  in  flux  the  anisotropies 
for  the  second  time  period  barely  exceeds  1.3  and  the 
temperature  ratio  never  rises  above  2.0. 

Figure  47  is  a  plot  of  the  logiof  versus  pitch  angle 
for  150  eV  electrons  from  different  days  and  hours  in  the 
detailed  analysis.  The  plot  is  organized  from  bottom  to  top  to 
show  a  progression  of  non-trapped  to  trapped  electron 
distributions.  Plot  a  and  b  show  the  pitch  angle  distribution 
for  periods  when  no  trapped  electron  distribution  is  present. 
Plot  c  and  d  show  the  beginnings  of  trapped  electron 
distributions.  Plots  e  and  f  show  the  narrowest  equatorially 
trapped  electron  distributions  observed  in  the  detailed 
analysis.  Note  that  as  the  trapped  distribution  gets  stronger 
the  peak  at  90°  gets  progressively  more  pronounced,  indicating 
that  the  electrons  are  equatorially  trapped.  This 
characteristic  suggests  that  the  data  may  be  fitted  with  a  bi- 
Maxwellian  distribution  function  as  was  done.  This  fit 
performed  well  and  it  was  quite  evident  it  did  not  fit  well 
indicating  that  the  distribution  was  not  trapped.  This 
criteria  allows  us  to  distinguish  between  (simply)  trapped 
distributions  (e.g.,  loss  c< nes) ,  versus  equatorially  trapped 
distributions . 
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A.  STATISTICAL  SURVEY 

The  probability  distributions  for  equatorially  trapped 
ions  and  electrons  show  clear  localization  for  regions  of  high 
occurrence  probability  (>=  50%) .  This  is  particularly  true  of 
electrons.  The  electron  distribution  has  a  very  localized  peak 
at  1000  “  1100  local  time,  L*=  6.0  -  6.5,  and  within  5  degrees 
of  the  magnetic  equator.  This  is  seen  regardless  of  the 
selection  criteria  used  in  the  surveys.  This  agrees  with  the 
results  found  by  Braccio.  The  survey  also  showed  the 
ec[uatorially  trapped  electrons  were  primarily  located  outside 
the  traditional  plasmasphere,  and  were  primarily  confined  to 
the  dawn  to  noon  sector.  The  peak  in  the  electron  distribution 
occurred  near  geosynchronous  orbit. 

The  results  for  trapped  ion  distributions,  also,  agree 
closely  with  Braccio 's  results  with  the  exception  of  his 
explanation  for  the  decrease  in  probability  between  1200  - 
1400.  He  suggested  that  the  decrease  in  probability  is  due  to 
the  average  energy  of  the  equatorially  trapped  ions  going 
above  the  150  eV  energy  he  surveyed,  in  this  local  time 
region.  Both  the  240  eV  and  442  eV  ions  surveyed  showed  the 
same  decrease  in  ion  probability  distribution  for  this  local 
time  region.  It  maybe  possible  that  Braccio 's  hypothesis  is 
correct.  Sagawa  et  al  (1987)  did  not  show  a  decrease  in 
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probability  for  the  region  from  1200  to  1500  local  time. 
However,  Sagava  et  al.  used  summary  results  from  approximately 
10  eV  to  1  keV,  whereas  Braccio's  survey  and  this  survey 
covered  50  to  442  eV.  To  accurately  determine  if  higher  energy 
ions  inhabit  this  region  of  decreased  probability  the 
additional  energy  channels  of  657,  885  and  1127  eV  ought  to  be 
surveyed.  The  ion  survey  reinforces  Braccios*  and  other 
previous  work  with  ISEE-1  and  DE-1,  with  peak  ion  flux  inside 
the  plasmapause. 

B.  DETAILED  ANALYSIS  CASE  STUDIES 

In  the  analysis  only  13  days  were  looked  at,  therefore  a 
statistical  analysis  isn't  appropriate.  However,  the  orbital 
coverage.  Figure  48,  of  the  days  analyzed  was  widespread 
enough  to  see  patterns  or  tendencies  emerge  in  the  data. 

Figure  49  is  a  plot  of  where  the  150  eV  data  from  the  13 
days  that  met  the  criteria  to  be  considered  a  trapped  electron 
distribution  occurred;  namely  a  flux  greater  than  5.0x10*  (cm^ 
s  sr)"^,  an  anisotropy  greater  than  2.0,  a  temperature  ratio 
greater  than  3.0  and  within  10®  of  the  magnetic  equator.  The 
data  is  plotted  in  the  x-y  plane  of  the  AMPTE  satellite  orbit. 
The  trapped  distributions  start  at  dawn  at  a  mean  L  value  of 
about  3.8  and  increase  in  orbital  density  and  L  to  a  maximum 
between  0900  and  1100  local  time  and  an  average  L  of  7.0.  The 
distributions  then  drop  off  rather  rapidly  with  only  smaller 
occurences  of  equatorial ly  trapped  distributions  in 
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mid-aftemoon  at  L  values  in  excess  of  7.0  and  at  dusk  when 
the  distributions  move  back  into  the  smaller  L  values  of  5.0 
to  6.5. 

Braccio  (1991)  in  his  statistical  survey  found  a  similar 
general  pattern  to  exist  for  trapped  electron  distributions  at 
the  150  eV  energy  range.  He  found  that  the  electron 
distribution  had  a  very  localized  peak  at  0900  local  time,  L 
=  6,  and  within  5°  of  the  magnetic  equator.  The  difference  in 
the  results  can  be  attributed  to  his  survey  considering  all 
days  of  operation  the  AMPTE  satellite  while  this  thesis  only 
considered  a  small  portion  of  the  operational  life. 

The  tendency  for  the  trapped  electron  distributions  to 
move  out  further  in  L  value  starting  at  dawn  until  late 
afternoon  and  then  move  back  in  near  dusk  can  be  explained  in 
terms  of  plasmapause  location.  Figure  50  is  taken  from 
Williams  et  al.  (1981).  It  shows  the  location  of  the 
plasmapause  as  a  function  of  local  time.  The  shape  of  this 
plot  mimics  the  behavior  of  the  trapped  electron  distributions 
in  the  case  studies.  Then  as  the  plasmapause  moves  out  during 
the  day  so  does  the  location  of  the  trapped  electrons. 

Braccio  (1991)  found  that  trapped  electrons  seem  to  be 
excluded  from  regions  of  high  trapped  ion  probabilities  and 
vice  versa.  He,  also,  found  that  the  trapped  electron 
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distributions  probabilities  drop  off  by  local  noon  at  the 
beginning  of  the  peak  of  the  trapped  ion  probabilities,  which 
he  found  to  peak  between  1400  and  1500  local  tine  at  L  =  8. 
Horowitz  et  ai.  (1981)  observed  that  trapped  ion  distributions 
tended  to  occur  along  the  plasnapause  boundary  in  the  dawn  to 
dusk  region  with  the  majority  just  inside  the  boundary. 
Therefore  Braccio  postulated  that  if  trapped  ions  occur  at  the 
plasnapause,  then  as  the  altitude  of  the  plasnapause  rises,  in 
the  course  of  a  day,  the  region  where  trapped  electron 
distributions  usually  occur  is  saturated  with  trapped  ions. 
This  ion  saturation  seems  to,  in  turn,  lead  to  a  disruption  in 
the  process  by  which  the  trapped  electrons  are  produced.  This 
approach  could  be  used  to  explain  the  gap  in  trapped  electron 
distributions  found  in  the  case  studies.  Figure  49,  between 
1200  and  1700  local  time  inside  of  L  »  7. 

Figure  51  is  a  plot  of  the  occurrence  of  340  eV  data  in 
the  13  days  that  met  the  criteria  to  be  considered  a  trapped 
electron  distribution;  namely  a  flux  greater  than  5.0x10®  (cm^ 
s  sr)‘^,  an  anisotropy  greater  than  4.4,  a  temperature  ratio 
greater  than  3.0  and  within  10°  of  the  magnetic  equator.  The 
data  is  plotted  in  the  x-y  plane  of  the  AMPTE  satellite  orbit. 
There  appears  to  be  no  major  grouping  of  trapped  electron 
distributions  at  this  energy  level.  This  is  probably  the 
result  of  the  higher  anisotropy  requirement  of  4.4.  The 
trapped  distributions  are  localized  to  0900  to  1100  local  time 
at  L  of  7.0  and  at  dusk  from  1700  to  1800  local  time  at  L=5.0. 
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A  comparison  of  the  two  different  anisotropies  used  to 
characterize  the  data  found  that  they  tracked  reasonedsly  well 
with  the  equator  anisotropy  being  on  average  1.075  times  the 
equator  anisotropy  in  the  dawn  to  noon  region  and  only  0.59 
times  the  equator  anisotropy  in  the  noon  to  post  dusk  region. 
The  difference  in  anisotropy  ratios  for  the  equator  anisotropy 
between  the  morning  and  the  afternoon  local  time  could  be  due 
to  a  lack  of  samples  taken  in  the  afternoon  region  that  were 
characterized  as  trapped. 

Over  28,000  samples  taken  in  13  orbits  are  summarized  in 
Figure  52,  the  plot  of  Tp„p/Tpaii  vs  Equator  Anisotropy  for  150 
eV  electrons.  This  figure  shows  a  different  technique  for 
defining  anisotropy  of  a  equatorially  trapped  electron 
distribution  and  suggests  that  there  is  a  monotonic 
relationship  between  a  temperature  ratio  greater  than  3  and  an 
anisotropy  of  2.0  for  150  eV  and  4.42  for  340  eV.  For  all  days 
the  threshold  of  3  for  a  temperature  ratio  begins  to  be 
exceeded  at  the  anisotropies  of  about  2.0  for  150  eV  and  about 
4.42  for  340  eV. 
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Figure  52.  Temperature  Ratio  verses  Equator  Anisotropy  -  All 
13  Detailed  Analysis  Case  Study  Days,  Energy  150  eV 
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V.  COHCLUSIOMS 


A.  STATISTICAL  ANALYSIS 

The  location  of  the  eguatorially  trapped  plasmas  are 
species  dependent.  The  ions  and  electrons  show  a  different 
local  time  dependence  in  the  location  of  their  occurrence 
probability  peak.  Electrons  show  a  uniform  high  probability 
distribution  centered  at  0900  local  time  and  an  L  value  of  6. 
The  fact  that  trapped  electrons  begin  to  be  seen  at  local  dawn 
lead  to  speculation  that  their  existence  is  dependent  on 
photoelectron  emission  form  the  Earth's  ionosphere.  The  shape 
of  this  distribution  is  basically  conical,  however,  it  drops 
off  more  rapidly  than  it  increases,  with  respect  to  local 
time. 

Ions,  meanwhile,  have  a  probability  distribution  that 
shows  a  strong  L  dependence  on  local  time.  The  high  trapped 
probability  region  begins  at  local  dawn,  for  L  approximately 
4,  and  rises  to  a  maximxim  at  between  1400  and  1500  local  time 
with  an  L  value  of  8.  The  distribution  then  drops  quickly  in 
altitude  as  local  time  increases.  Additionally,  there  is  a 
region  of  decreased  probability  in  the  afternoon  sector  for 
this  data  that  was  speculated  by  Braccio  (1991)  to  be 
inhabited  by  higher  energy  ions,  but  still  is  present  for  the 
higher  energies  surveyed. 

The  overall  shape  of  the  high  trapped  ion  probability 
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region  mirrors  that  of  the  location  of  the  plasmapause.  This 
suggests  that  the  trapped  ion  distribution  is  linked  to  the 
plasmapause.  This  would  confirm  Horwitz  et  al.  's  (1981) 
observations  that  'pancake'  distributions  often  occur  in  the 
vicinity  of  the  plasmapause. 

Additionally,  it  has  been  observed  that  trapped  electrons 
seem  to  be  excluded  from  regions  of  high  trapped  ion 
probabilities  and  vice  versa.  If  the  trapped  ions  actually  do 
occur  at  the  plasmapause,  then  as  the  altitude  of  the 
plasmapause  rises,  in  the  course  of  a  day,  this  exclusion 
would  effectively  create  a  barrier  that  restricts  the  trapped 
electron  distribution  to  the  dawn  to  noon  sector. 

B.  DETAILED  ANALYSIS  CASE  STUDIES 

The  results  of  the  detailed  analysis  case  studies  in 
general  agree  with  the  above  results  found  by  Braccio  (1991) . 
Trapped  electron  distributions  are  concentrated  in  the  dawn  to 
noon  ragion  with  a  peak  around  1000  at  an  L  value  of  7.  What 
trapped  distributions  occur  in  the  afternoon  are  generally  at 
higher  L  values  until  late  afternoon  when  they  move  back  into 
L  values  around  5  to  6  at  1800.  This  follows  with  the  results 
found  by  Braccio  (1991)  who  found  that  large  trapped  electron 
distributions  are  exclude  from  areas  of  high  trapped  ion 
distributions  that  peak  in  the  1400  to  1500  local  time  region. 
This  also  suggests  a  connection  between  the  plasmapause  and 
trapped  electron  distributions. 
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The  threshold  flux  level  used  in  the  statistical  survey 
was  too  low  for  150  eV  electrons.  A  higher  flux  level  of  5x10' 
(cm^  s  sr)'^  could  be  used  for  this  energy  level.  The  use  of  an 
anisotropy  of  1.5  and  a  temperature  ratio  of  2.0  produced 
results  more  in  keeping  with  the  results  of  the  statistical 
survey.  The  threshold  flux  level  for  340  eV  appeared  to  be 
correct . 

The  use  of  an  average  anisotropy  at  a  latitude  within  6° 
of  the  magnetic  equator  vice  a  calculated  equator  anisotropy 
in  defining  trapped  electron  distributions  is  generally  pretty 
accurate.  The  average  difference  is  only  7.5%  between  the  two. 

It  has  been  shown  that  equatorially  trapped  plasmas  can  be 
described  by  a  bi-Maxwell ian  distribution  function.  With  such 
a  description  the  values  of  the  perpendicular  and  parallel 
temperatures  can  be  calculated  and  their  ratio  used  along  with 
an  anisotropy  ratio  between  the  flux  in  the  80®  to  90°  pitch 
angle  range  with  the  flux  in  the  60°  to  70°  pitch  angle  range 
to  define  a  electron  distribution  as  equatorially  trapped  or 
not.  It  was  found  that  a  temperature  ratio  of  three  along  with 
anisotropies  of  2.0  for  150  eV  electrons  and  4.4  for  340  eV 
electrons  could  be  used  to  define  a  trapped  electron 
distribution  along  with  the  required  flux  and  magnetic 
latitude  requirements. 

However,  the  340  eV  anisotropy  of  4.4  produce  less 
reliable  results  than  the  150  eV  anisotropy  with  the 
theoretical  anisotropy  and  temperature  ratios  for  340  eV  being 
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too  high.  In  this  study  Tp^^p  is  calculated  from  the  bi- 
Maxwellian  distribution  function  and  this  results  in  the  same 
Tp,rp  for  both  150  eV  and  340  eV.  This  was  the  only  method 
available.  However,  the  method  produces  a  Tp^^p  that  is  an 
underestimate  of  the  perpendicular  temperature  for  340  eV 
particles.  This  suggests  that  the  trapped  electron 
distributions  at  this  energy  are  not  true  bi-Maxwellians,  but 
are  more  like  a  kappa  function. 


APPENDIX  (DETAILED  ANALYSIS  SPECTROGRAMS) 
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Figure  53.  Pitch  Angle  verses  Time  Spectrogram  -  Day  84/266, 
Energy  150  eV 


112 


AMPTE  84283  150  eV  Electrons 


OT  4  0  4fi  S.0  &5  BO 

LT  52  77  SB  04  09 


117  2.9  &0  4.7 


in 

•  0 

•  9 

00 

e  9 

10  0 

LT  10^ 

tl.t 

11-9 

11.4 

11.9 

L 

7.7 

«.o 

a.3 

0.9 

0.7 

m  10.0 

10.9 

11.0 

11.9 

U  11.9 

11J 

11.0 

12.1 

t  a.7 

08 

0.9 

09 

X  -^0 

.  .  .  . 

_ _  -9  *  .  . 

0  4 

j 

c’~. 

$  M-  .  . 

Figure  55.  Pitch  Angle  verses  Time  Spectrogram  -  Day  84/296, 
Energy  150  eV 
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Figure  57.  Pitch  Angle  verses  Time  Spectrogram  -  Day  85/002, 
Energy  150  eV 
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Figure  58.  Pitch  Angle  verses  Time  Spectrogram  -  Day  85/039 
Energy  150  eV 
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Figure  59.  Pitch  Angle  verses  Time  Spectrogram 
Energy  150  eV 
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Figure  60.  Pitch  Angle  verses  Time  Spectrogram  -  Day  85/238 
Energy  150  eV 
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Figure  63 .  Pitch  Angle  verses  Time  Spectrogram 
Energy  150  eV 
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Figure  68.  Pitch  Angle  verses  Time  Spectrogram  -  Day  84/315, 
Energy  340  eV 
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Figure  70.  Pitch  Angle  verses  Time  Spectrogram 
Energy  340  eV 
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AMPTE  85264  340  eV  Electrons 


Figure  75.  Pitch  Angle  verses  Time  Spectrogrcim  -  Day  85/301 
Energy  340  eV 
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Figure  76.  Pitch  Angle  verses  Time  Spectrogram  -  Day  85/320 
Energy  340  eV 
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Figure  77.  Pitch  Angle  verses  Time  Spectrogram  -  Day  85/337 
Energy  340  eV 
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